use "BJPS_WVS_Data.dta", clear




*** FIGURE A9
*** TRUST OTHER REL  
reg TrustOtherRel i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store A

*** RELIGIOUS FAITH IMPORTANT IN CHILD  
reg ChildFaithNotImportant i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store B

*** CONFIDENCE IN CHURCHES  
 
reg NotConfidentInChurches i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store C

*** ONLY ACCEPTABLE RELIGION  
reg OtherReligionsAcceptable i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 [pweight=S017] 
margins, dydx(Friday) over(Weekly) post
estimates store D

coefplot A B C D, xline(0) scheme(s1mono) legend(position(6) cols(1)) graphregion(color(white)) bgcolor(white) plotlabels("Trust other religions" "Child faith not important" "Not confident in religious institutions" "Other religions acceptable") levels(95 90) xtitle(Marginal Effect of Friday Interview)


*** FIGURE A10
**** ONLY NONMUSLIM MAJORITY COUNTRIES
*** TRUST OTHER REL 
reg TrustOtherRel i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim < 0.5 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store A

*** RELIGIOUS FAITH IMPORTANT IN CHILD   
reg ChildFaithNotImportant i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim < 0.5 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store B

*** CONFIDENCE IN CHURCHES 
reg NotConfidentInChurches i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim < 0.5 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store C

*** ONLY ACCEPTABLE RELIGION 
reg OtherReligionsAcceptable i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim < 0.5 [pweight=S017] 
margins, dydx(Friday) over(Weekly) post
estimates store D

coefplot A B C D, xline(0) scheme(s1mono) legend(position(6) cols(1)) graphregion(color(white)) bgcolor(white) plotlabels("Trust other religions" "Child faith not important" "Not confident in religious institutions" "Other religions acceptable") levels(95 90) xtitle(Marginal Effect of Friday Interview) saving(nonmuslim_countries.gph, replace) title(Non-Muslim-Majority Countries) yscale(alt)


	**** MUSLIM-MAJORITY COUNTRIES 
*** TRUST OTHER REL
reg TrustOtherRel i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim > 0.5 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store A2

*** RELIGIOUS FAITH IMPORTANT IN CHILD  
reg ChildFaithNotImportant i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed  i.S002 if PctMuslim > 0.5 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store B2

*** CONFIDENCE IN CHURCHES 
reg NotConfidentInChurches i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim > 0.5 [pweight=S017]
margins, dydx(Friday) over(Weekly) post
estimates store C2

*** ONLY ACCEPTABLE RELIGION 
reg OtherReligionsAcceptable i.Friday##i.Weekly i.S003 X025R X003 X001 PrayDaily Unemployed i.S002 if PctMuslim > 0.5 [pweight=S017] 
margins, dydx(Friday) over(Weekly) post
estimates store D2

coefplot A2 B2 C2 D2, bylabel(Muslim Countries) || A B C D, bylabel(Non-Muslim Countries) ||, xline(0) scheme(s1mono) graphregion(color(white)) bgcolor(white) plotlabels("Trust other religions" "Child faith not important" "Not confident in religious institutions" "Other religions acceptable") levels(95 90) xtitle(Marginal Effect of Friday Interview)  legend(cols(1))

 